A note on the critical nuclear charge that binds iV electrons 
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The critical charge, Z c that binds TV electrons in an atom is known to satisfy the relation, Z c — 
TV — 1 — 8. Numerical calculations show that 5 exhibits a correlation with the electron affinity, E a , 
of the neutral atom with TV — 1 electrons, and that 5 < 0.48 for TV < 90. We re-examined Z c by 
means of a very simple model that treats the interaction of the outermost electron with the core 
electrons through a zero-range force theory, whereas the residual Coulomb interaction is included in 
first-order perturbation theory. The relation, S = 0.836 V E a , coming from our model, is fairly well 
observed by the computed S and the experimental E a . 

PACS numbers: 32.30.-r, 32.10.Hq, 31.15.-p 



I. INTRODUCTION 

Qualitative understanding of phenomena come from 
very simple models that capture their essence. In atoms, 
perhaps the most successful simple models are the fa- 
mous Thomas-Fermi (TF) Theory^l, and the large-D 
expansion^, leading, the latter, to a qualitative picture 
similar to the one devised by Lewis long ago3J 

The theory by Thomas and Fermi is a mean-field one 
in which the scaling properties of the kinetic energy and 
Coulomb potentials are apparent. This fact leads to a 
universal relation for the energy of the atom with TV 
electrons and charge Z®E TF (N,Z) = N 7 / 3 f(N/Z), or 
for the ionization potentiaP, I P (N,Z) = N 4 / 3 g(N/Z), 
where the universal functions / and g depend only on 
the combination N/Z. 

For a fixed ratio N/Z, the ionization potential (or the 
electron affinity, E a ) should behave as E a ~ TV 4 / 3 , when 
N increases. This dependence holds extremely well for 
cations (TV < Z), even for small TV values, as can be seen 
in the upper panel of Fig. [TJ but breaks down for neutral 
atoms (TV = Z), where I p or E a , apart from fluctuations 
due to shell effects, take roughly a constant value (Fig. 
[I] lower panel). 

The reason for such a breaking down of TF predictions 
for neutral atoms, in certain situations, is the proximity 
of the anionic instability threshold, that is, the fact that 
atoms become unstable when Z < Z c < TV — 1. 

Rigorous bounds for the critical nuclear charge, Z c , 
that binds TV electrons have proven to be of little 
use^. N umer ical computations for atoms with up to 89 
electronJ^^, on the other hand, lead to: 
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FIG. 1: (Color online) Upper panel: Electron affinities of 
cations with TV < 20 and fixed N/Z ratios. Data are taken 
from the NIST database® The line 0.4 N 4/3 is drawn as a 
reference for the slope. Lower panel: Electron affinities of 
neutral atomS 



Z r = TV - 1 



(1) 



where the magnitude S < 0.48 shows a correlation with 
the electron affinity of the neutral atom with TV — 1 elec- 
trons. 

In the present paper, we use a very simple model, si- 
milar in spirit to the one used in Refs. |10j . from which 
we obtain (in atomic units): 



S = 0.836a 



(2) 



Relation ^ is fairly well observed by the calculated 
S when the experimental E a of the neutral atoms with 
TV — 1 electrons are introduced, as shown in Fig. [2] 



II. THE MODEL 

When Z « TV — 1, we may use a very simple model to 
get reasonable results. The interaction of the outermost 
electron with the neutral core is described by a square 
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FIG. 2: (Color online) Computed values of 8 versus 
the square root of the experimental electron affinity for the 
neutral atom with N — 1 electrons-^. The central line is our 
prediction, Eq. Q 



well potential of depth V and radius R. V and R enter 
in a combination such that the binding energy of the 
electron is E a . 

At this level, the parameter R is not determined, as 
in a theory with zero-range forces^, and the details of 
the potential for r < R are not important. The only 
magnitude that matters is E a , which determines the wave 
function in the outer region: 



'2k 
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(3) 



where k = \j2E a . This function can be used to pertur- 
batively compute the contribution to the energy coming 
from the Coulomb interaction with the almost neutral 
core: 



f°° dr 

AE = (N - 1 - Z) 2k I — exp(~2nr) 
Jr r 



{N-l-Z) \jE a F(\jE a R), 



(4) 



where the explicit expression F{x) = 2v / 2(ln l/(2y/2x) — 
0.5772) +8x comes from the assumption that \/k should 
be greater than the radius R. 

The energy of the anion with N electrons and charge 
Z m N - 1 is, thus: 



E = -E a + (N-l-Z) ^E a F{^E a R). (5) 

We may use the freedom left in the definition of the 
parameter R in order to minimize the energy. It may be 
seen that F takes a minimum value, 1.196, when R = 
0.35/ jE~ a (Fig. U upper panel). 

When E a is small enough, one may expect that the 
slope of the energy curve at Z = N — 1, i. e. —1.2 y/Ea, 
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FIG. 3: (Color online) Upper panel: The function F entering 
Eq. Q. It is computed in lowest order, up to terms linear in 
\SEaR- Lower panel: Comparison between the values taken 
by the parameter R = 0.35/\AEa, the covalent radius, and 
the radius used by Kais and co-workers^ for the neutral atom 
with N — 1 electrons. 



is practically the same as the slope computed at Z c . The 
linear behaviour at Z c was demonstrated in Ref. [12] . 
Thus, we may use this slope to estimate Z c , leading to 
the expression for S given in Eq. ^ of section [I] 



III. DISCUSSION 

We put together in Fig. [2] the computed values of 
S and the experimental EJ^. Most of the points are 
grouped around the line given by Eq. ((2). There are. 
however, a few exceptions that we want to comment. 
They may correspond to inconsistencies in the computed 
5 or the measured E a . 

The point with 7Y = 86 is the only one for which E a 
is not measured, but computed^. The reported electron 
affinity of At seems to be underestimated. 

The points with N — 60 and 70 are related to the E a 
of Pr and Tm, respectively. Both measurements were 
performed by the same group using the same experimen- 
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tal seturj^. The authors stress that the reported val- 
ues, around 1 eV, are much higher than the theoretical 
estimates^, around 0.1 eV. We notice that a value of 0.1 
eV is more consistent with our prediction. 

On the other hand, the E a of Fe and Ti enter the points 
in the Figure labelled as N = 27 and 23. We notice that 
there is also a small difference between measured-^ and 
theoretical values. The latter are closer to our curve. 

Let us stress that the radius R = 0.35 / \/Ea takes val- 
ues in between the covalent radius and the radius used 
by Kais and co-worker for the neutral atom with N — 1 
electrons, as shown in Fig. [3j lower panel. The excep- 
tions are the systems with very low E a , for which R be- 
comes higher than the covalent radius. The magnitude R 
is simply a parameter in our model, not directly related 
to the real potential that experiences the electron. It is 
remarkable that the theory, in the kR — > limit, assigns 
to R "reasonable" values. 



IV. CONCLUSIONS 

We use a kind of "zero-range force theory" in order 
to compute, for the atom with N electrons and charge 



Z, the slope of the curve E vs Z at the point where 
Z = N — 1. The obtained slope, 1.2 \[W a , allows us to 
estimate the critical charge for binding N electrons, Z c = 
N — 1— 0.836\/E^. Comparison between the computed Z c 
and measured E a shows that this simple relation holds. 
Some points lying far from the curve are characterized by 
discrepancies between the experimental and theoretical 
values for E a . 

The case E a = 0, Z c = N — 1 is also contained in 
our curve. There is an additional property of the energy 
curve in this situation. Due to the fact that it is precisely 
the residual Coulomb interaction which binds the elec- 
tron, the energy should behave as: E ~ —(Z — N + l) 2 , 
for Z close enough to iV — 1. We did not find computed 
data to check this behaviour. 
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